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Abstract 


Large- scale buoyant low stretch stagnation point diffusion flames 
over a solid fuel (polymethylmethacrylate) were studied for a range of 
aerodynamic stretch rates of 2-12 sec' 1 which are of the same order as 
spacecraft ventilation-induced stretch in a microgravity environment. An 
extensive layer of polymer material above the glass transition temperature 
is observed. Unique phenomena associated with this extensive glass layer 
included substantial swelling of the burning surface, in-depth hubhU 
formation, and migration and/or elongation of the bubbles normal Uj the- hui 
surface. The bubble layer acted to insulate the polymer surface by reducing 
the effective conductivity of the solid. The reduced in-depth conduction 
stabilized the flame for longer than expected from theory neglecting the 
bubble layer. While buoyancy acts to move the bubbles deeper into the 
molten polymer, thermocapillary forces and surface regression both act to 
bring the bubbles to the burning surface. Bubble layers may thus be very 
important in low gravity (low stretch) burning of materials. As bubbles 
reached the burning surface, monomer fuel vapors jetted from the surface, 
enhancing burning by entraining ambient air flow. Popping of these bubbles 
at the surface can expel burning droplets of the molten material, which may 
increase the fire propagation hazards at low stretch rates. 


2 



simulate their combustion behavior in low gravity. The first space-based 
experiment on burning solids was conducted aboard Skylab in 1974 iU . Tin- 
most dramatic test in this series involved a piece of nylon which burned for 
almost 11 minutes in a quiescent 65% 02, 5.2 psia atmosphere before tho 
astronauts extinguished it via evacuation of the chamber. During the te.-t. 
the sample melted and coalesced into a molten globule. Boiling of the high 
viscosity molten sample was noted, and it burned with an agitated, 
pulsating flame. 

Although imaging was poor, the sooty flame pulsations appeared to 
be due to local nylon fuel vapor jetting of the boiling melt. Paper samples, 
which were also tested during the Skylab flight, extinguished after flame 
spread was completed but before the sample was consumed. The 
extinguishment was attributed to the accumulation of combustion products 
and insufficient oxygen to the flame ll *. The contrast between the nylon and 
paper flames is revealing: a flame in a quiescent environment will self- 



extinguish duo to accumulation of combustion products , but self-induced flow- 
disturbances from the boiling/ vapor jetting are sufficient to sustain a flame. 

Other low gravity testing with nylon velcro 1 " 1 and 
polymethylmethacrylate (PMMA) spheres 1 * 1 also revealed vigorous perturbed 
burning of the molten materials due to fuel vapor jetting. Similar to the 
paper results above, recent testing in a Space Shuttle experiment l; 
demonstrated self-extinguishment after more than 9 minutes oi burning ui .1 
small metal-housed PMMA slab in a quiescent, 50% 02. 1 atm microgravity 
environment. Prior to flame extinction in that Space Shuttle test, one 
instance of apparent vapor jetting is noted, and the flame stabilized for an 
extended period of time (flow decay time scales) due to that event. However, 
due to the small sample size and high heat loss to the housing, the solid 
material did not melt during the test. As such, the atmosphere was 
essentially quiescent, with the accumulation of combustion products reducing 
the local oxygen concentration, causing extinction. Thus it appears that 
materials that melt and have significant vapor jetting as they burn pose a 
more significant spacecraft fire hazard due to the self-induced flow they 
generate. 
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The characteristics of molten layers of burning materials are nm 
well understood even in normal gravity. Visual observations of thin bubble 
layers have been reported for normal gravity PMMA burning, but with very 
little in the way of supporting analysis 1 ’ , ’ 1 . The importance of the bubble 
layers on heat transfer processes such as heat conduction and 
absorption/scattering of incident external radiation has not been examined 
Thick molten layers with copious in-depth fuel vapor production may play a 
significant role in the viability of flames in microgravity. This study is in-.- 
first to look at the characteristics and influences of the molten and bubble 
layers formed in PMMA burning at low stretch rates (low heating rates). The 
low stretch rates are of the same order as the ventilation-induced stretch 
rates in a microgravity environment such as on a spacecraft. In this article, 
we report observations of the characteristics of the molten and bubble 
layers, and provide some scaling of the bubble layer’s importance in heat 
and mass transport during the burning process. Details on the gas-phase 
studies are reported in [8,9], This paper focuses on the solid phase since 
important polymer degradation characteristics are not nearly as well 
understood as the gas-phase processes. 
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Polymethylmethacrylate (PMMA) Characteristics 


Commercial cast PMMA is an amorphous, uncrosslinked polymer of 
methyl methacrylate, the monomer. The monomer units are oriented in the 
chain in an atactic to syndiotactic configuration so that the large ester group 
is either randomly oriented or alternates sides in adjacent monomers 1 . 
Chain lengths are very long, typically on the order of thousands of monomer 
units. PMMA was the first polymer that was fuiind tu break up upon 
pyrolysis into almost entirely monomer units by a stepwise unzipping process 
[U1 so that the fuel vapor is predominantly methyl methacrylate (MALA). 

This is believed to happen by initially random scissions of chains. Due to 
steric hindrance from the methyl and ester groups on the quaternary carbon 
from the main chain, hydrogen transfer doesn't occur and the chain break 
point ends with two free radicals. Each end of the chain then forms a double 
bond within the monomer backbone, which detaches the monomer from the 
chain and creates a new free radical end. As such, the free radicals unzip like 
a broken string of pearls to yield double-bonded monomer units. Table 1 
contains PMMA material characteristics information. 
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Various vendors make cast PMMA, and there are differences between 


vendors, primarily in the average molecular weights of the resultant polymer 
chains. The commercial casting process is a hatch polymerization of a syrup 
of monomer/polymer in the presence of small quantities of the initiator 
benzoyl peroxide. Since the process is diffusion controlled, the chains are 
“tangled" and show no directionality. In contrast, extruded PMMA has a 
distinct directionality of the chains due to the continuous flow 
polymerization. Because of their "untangled" orientation, the extruded 
samples have a much greater tendency to drip while burning since the 
material has a low glass transition temperature relative to the pyrolysis 
temperature. 

Due to cast PMMA's desirable burning characteristics ( it degrades to 
monomer, its degradation kinetics are relatively well understood, and it does 
not char or drip), it has become a classic non-charring fuel used in solid-fuel 
combustion studies. Much of what we currently understand about the 
combustion of thermally-thick materials is derived from testing with cast 
PMMA. [5 12,131 . Most studies treat the material as a solid that pyrolyzes to a 
gas only at the burning surface via a zeroth order Arrhenius kinetic model 
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Initial models with subsurface degradation assumed instantaneous transport 
of the gaseous products to the surface* 11 *. Solid-phase properties are tvpicallv 
assumed to be constant throughout the condensed phase (density, thermal 
conductivity, thermal diffusivity. heat capacity). 

However, these assumptions may be misleading, especially for low 
stretch flames close to an extinction limit. A two-phase layer at the burning 
surface of PMAIA has been noted repeatedly . The soiid phase uu.- 
observed to consist of the solid, a molten polymer layer, and then a surface 
two-phase (bubble) layer, across which a temperature change of 
approximately 300 K is realized. Subsurface thermocouples were noted to 
become unreliable as they approached the burning surface, due to the effects 
of a 'liquid-like layer' adjacent to the surface of the PMMA. Also noted in 
these studies was variability in the extinction limits during testing, due to 
"random disturbances", which was interpreted as being due to small bubbles 
in the two phase layer that rupture and release a vapor jet when they reach 
the surface. 


Vovelle et al. * 1SI found in their thermal degradation (not combustion) 
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experiments with PMMA samples under steady external radiant heat flux 
that the sample size must he at least 1.7 to 1.<S cm thick to obtain a 
consistent maximum rate of mass loss, although mass loss varied in time 
even for the thickest samples. Thinner samples provided increased heat loss 
from the backside of the specimens, which lowered the maximum rate of 
mass loss achieved. The solid phase was always unsteady over long times 
(on the order of 1 hour)., as measured by the mass loss rate. The subsurface 
region (a two-phase layer) contributed to the rate of gasification. 

Chaiken et al 1191 found for surface temperatures below 714 K (typical 
PMMA burning surface temperatures are nearer 650K), that the polymer 
surface is not saturated with monomer, and its formation and diffusion to 
the surface are the rate-controlling processes of gasification. Thus most 
combustion situations with PMMA have an unsaturated polymer surface, 
and the in-depth monomer formation and diffusion are the rate-controlling 
processes in the pyrolysis process. Heat absorbed from the flame in this 
layer, in the form of conduction and radiation, directly influences pyrolysis. 

It is important to understand the behavior of the near surface layer of 
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burning PMMA, since it plays such a significant role* in the vaporization 
process. However, no detailed experimental results have been reported of 
this layer. This paper presents some data on the near surface layer of PMMA 
burning under low stretch conditions. 

Experiment 


A noimal giavit\ low stretch flame staoili/cd over a sulid fuel iici.-? oeeii 
studied through the use of large radii samples in a stagnation point 
geometry. For purely buoyant flow, a diffusion flame can be established at 
the bottom stagnation point of a blunt body. In this case the stretch rate, 
a * (g/R)*, where g is gravity level and R is the radius of curvature of the 
body. 


The cylindrical segments used as fuel samples in these experiments 
were made from cast PMMA stock sheets. Sections roughly 30 cm (12”) wide 
by 30 cm (12") long by 2.4 cm (0.944 -1 ) thick were cut from the flat stock 
(oversized slightly to allow for 2% shrinkage). These samples were heated 
under infrared lamps for 90 minutes to a forming temperature near 160°C. 
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At that temperature, which is well above the "lass transition temperature (T 
: 105"C), the samples arc rubbery-soft and pliable 1 . The hot samples were 
placed within a wood mold with felt covers and press-shaped in the desired 
radius of curvature for the cylindrical sections. They were insulated within 
the mold and left to cool slowly to allow for stress relaxation during tin- 
cooling process. Samples were then cut to size and sample curvature was 
measured manually to confirm the correct radii. Samples made this way 
varied from 20-200 cm radius. The 5 cm radius samples were cut from a 
cast solid cylinder into half-cylinders, due to the difficulty in forming such 
tight radii with the above technique. 

The samples were ignited uniformly over their lower surface using a 
gas-jet flame array, and a flame was established below the sample. Ignition 
times valued from 1-15 minutes, with larger radii needing more time to 
ignite. The ignition flame array provides more heat flux than the stable 
flame in order to heat the cold sample and initiate surface burning. The 
influence of the ignition on measured swelling, temperature distributions, 
and bubble formation is negligible, since the swelling, temperature 
distributions, and bubble layers continue to develop after the igniter is 
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deactivated, as will he shown in the Sample Swelling and Regression section. 
Figure 1 shows a conceptual schematic of the burning sample, as seen along 
its axis. The flame obtained after an initial solid-phase heatup period is a 
one-dimensional, uniform flame, with a pseudo-constant surface-to-flame 
standoff distance in the stagnation zone of the cylinder. The solid-phase is 
inherently transient (over long times) because of surface regression of the 
sample’s finite thickness. The burning material is cooled on the top surface 
with an ice bath, and the solid-phase temperature increases smoothie 
(though not linearly as will be discussed later) through the sample to a 
pyrolysis temperature of approximately 630K at the regressing surface. The 
surface is somewhat uneven during the burn, due to local dripping in some 
instances, and due to vapor bubble rupture at the surface in others. A 
bubble layer is formed near the burning surface, and deeper within the 
PMMA sample a refractive index change surface is noted in the post- 
combustion samples. Different radii samples (different stretch rates) were 
examined to determine trends in the these features. 

Thermocouples, embedded at different depths (four 0.0127 cm OD 
bare wire type K, with leads normal to the burning surface) were used to 
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record solid temperatures prior to and during selected tests. The 
thermocouple readings were sampled at 1/(5 Hz, signal conditioned and cold- 
junction compensated, with an uncorreeted accuracy of 20 0. An infravd 
radiometer with a Sehmidt-Boelter thermopile sensor (0.0-15 microns) was 
used to record the total radiant loss from the flame during the experiment. 


Burning rates, gas and solid-phase temperatures, and radiation loss 
were measured during the experiment so that detailed surface energy 
balance calculations could be made. The full details of these experiments 
are reported in Olson [8 91 . The samples were allowed to burn from 5 minutes 
to 3 hours, depending on their burning rate (radius of curvature, or 
equivalently the stretch rate). At the end of each test, the flame was quickly 
extinguished with a manually-applied water mist/spray, which did not 
disturb the surface significantly. The burned PMMA was allowed to cool 
completely before being moved. Temperature readings in the sample were 
below the glass temperature within a few minutes of extinguishment. 

Surface and Bubble Layer characteristics 
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Post-burn samples were cut normal to the cylinder axis and polished 
to examine the character of the solid phase preserved from the burn. 
Photographs were taken of the front surface, hack surface, and huhlde Inver 
Measurements were made from the sample itself and these photographs. A 
summary of these measurements is shown in Table 2. 

The burned surface characteristics are shown in Figure 2 for o. 20. 50. 
and 100 cm radius-of-curvature samples. In Fig. 2a, the burned surface wa.- 
black with very uniform small craters approximately 1/3 mm in diameter. 
The black color is due to soot deposition onto the molten surface via gas- 
phase thermophoretic forces on the soot particles. In Fig. 2b, the surface of 
the faintly caramel-colored translucent burned sample was less uniform than 
Fig. 2a, with larger, more varied craters (0.5-2 mm in diameter). This 
surface was mottled with black globules of re-solidified PMMA, where the 
blackening agent is likely to be soot and caramel-colored tar-like pyrolysis 
products from the flame. These globules are frequently hollow and range 
from 0.5-1 cm in characteristic dimension. In Fig. 2c, the burned surface is 
more glassy in appearance with a darker caramel color on the surface. The 
size distribution of the surface craters is much larger than in Fig 2b, from 
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sub-mm to larger than a mm in a few instances, awl protruding from the 
substrate up to 3 mm. The drip globules are still black. As the radius 
increases in Fig. 2d. these trends are amplified, where the surface has ■> 
deep caramel color and the crater distribution larger. 


The bubble layer characteristics are shown in Figure 3 for 5, 20. 50. 
and 100 cm radius of curvature samples. In Fig. 3a. the bubble layer 
extended approximately 3.5 mm below the surface. The bubbles adjacent to 
the surface were typically small, glassy, and closely spaced, in agreement 
with the crater sizes noted in Fig 2a. The bubbles were frequently attached 
to the burned surface via a glassy neck. Deeper in the sample, larger near- 
spherical bubbles were predominant, with sizes up to 2 mm in diameter. 

Some were a transparent glassy color, while others were opaque white. In 
Fig 3b, the bubble layer is approximately 4.5 mm thick and is characterized 
by an extremely uniform dense distribution of elongated bubble tubes. The 
amber color of the uneven surface is not carried in-depth, and all of the 
bubbles have a transparent glassy appearance. The deepest bubbles are 
typically as large or larger than near surface bubbles (up to 1 mm), and have 
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a more spherical shape than those near the surface. Near surface bubbles 
are elongated and frequently joined to other bubble tubes. The diameter of 
these tubes is approximately 0.4 mm. with sub-mm sparine between tuhr~ 

In Fig. 3c, the thickness of the bubble layer is 5 mm. Unlike Fig 3b 
which showed a highly uniform and dense bubble layer. Fig 3c has a less 
dense surface bubble layer, with generally larger (0.5-1 mm diameter) and 
less elongated glassy bubbles. The caramel color also seems to permeate [ jr.- 
near surface layer. Deeper in the sample, there are a substantial number of 
large elongated crumpled opaque bubbles. They are typically tethered to 
the burned surface with a glassy neck. There are also a few spherical 
bubbles suspended within the sample, which are typically at least partially 
glassy with only sections appearing to be wrinkled and frosted. In Fig. 3d, 
the bubble layer is deeper (0.7 cm), with larger crumpled bubbles. For all the 
bubble layers, the orientation of the tubes and elongated bubbles is 
orthogonal to the surface tangent rather than strictly along the gravity 
vector, indicating that surface tension, the temperature gradient, and the 
viscous nature of the solid are also important parameters in the bubble layer 
development. 
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The source of the hubbies is believed to be in-depth random scission of 
the polymer chains and vaporization of the released monomer. This 
temperature dependent process, which starts at temperatures ns low ns 
20CTC, is most likely to occur near the burning surface (3G0 C). but it can 
occur in depth as well, especially if conduction heat transfer is augmented In- 
in depth absorption of the gas-phase radiation |JI)| . The glassy bubbles are 
simply residual MMA vapor (b.p. 101°C @ 1 atm.) pockets in the re-solidified 
amorphous polymer. The opaque, wrinkled bubbles are the same ALMA vapor 
pockets where the opaque wrinkling could be due to local fracturing of the 
cooling/solidifying molten polymer-vapor interface as the pressure within the 
vapor bubble decreases after the experiment due to cooling/condensation of 
the MMA vapor. 

Deeper yet within each solid (=> 1 cm), a transparent surface is 
visible which parallels the burned surface in the post-burn samples. This 
optical anisotropy (shown schematically in Fig. 1) is not crazing (i.e. a fine 
crack), but is believed to be related to the deepest penetration of the glass 
transition surface. The originally stress annealed sample undergoes self- 
weighting stress of the molten polymer at temperatures above the glass 
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transition, and this stress is frozen into the sample during the sudden 
quench at the end of the test. This change in stress history results in a 
change in the stress optical coefficient of the pnlvmmv 

The Importance of Solid-Phase Peclet Numbers 


To further understand the bubble layer within the solid, it is helpful t 
understand the temperature distribution within the material, f ur sanipio 
with the radius of curvature R much greater than the sample thickness L 
(where L is an intermediate thickness that is neither thermally thick or 
thermally thin), a one dimensional steady state conduction-convection 
equation with appropriate boundary conditions can be analytically solved as 
a function of depth and Peclet number, where the Peclet number is a ratio of 
convection (surface regression) to conduction, Pe— L V„ / a. The solution is 

rcLji 


0(P) = exp PfK - exp 
1 - exp' p< ' 


(Eqn. 1) 


where @(Y)=(T-TJ/(T S -T V ) is the non-dimensional temperature, and Y=j7L is 
the non-dimensional depth from the burning surface. L is the thickness of 
the sample, y is the distance orthogonal to the burning surface, a is the solid 
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phase thermal diffusivity (1.2x10 ! cnr/s), T„ is the surface temperature, and 
T is ambient temperature. 


As shown in Figure 4. a linear temperature profile is obtained for tin- 
limiting case of Pe=0, and the profile curvature increases with increasing Pe 
(increasing importance of regression). The thickness of the glassy layer (6 
line in Fig. 4), which is the layer of material heated above the glass 
transition, increases as Peclet number decreases. For Pe=0, the glassy layer 
is over 75% of the thickness of the sample. For Pe=10, on the other hand, 
only 15% of the sample is glassy. Peclet numbers in these experiments 
varied from 3.6 to 0.15, with the majority less than unity. Clearly, the glassy 
(i.e. molten) layer becomes dominant at low regression rates. 

Also shown in Figure 4 are the experimental measurements (for the 
three samples where the model’s R»L assumption is valid, i.e. 20, 50, and 
100 cm radius data from Table 2) of the bubble layer thickness and the 
position of the index of refraction surface, using the measured Peclet number 
(see Table 2) to estimate the non-dimensional steady-state temperature 
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using Eqn. 1. Although the steady state approximation is not rigorously 
valid due to regression of the finite thickness sample, it is a reasonable 
approximation if the final thickness of the sample is used to nemvilv". v v 
The bubble layer starts at non-dimensional temperatures as low as -().(*» 
(225"C), which agrees with data indicating decomposition becomes 
appreciable by about 200°C |l ". The index of refraction location is very close 
to the glass transition location, as previously suggested. 

Due to the strong curvature in the temperature profiles at high Peclet 
numbers, the temperature gradients at the regressing surface are quite 
steep. In contrast, the temperature gradients at the back surface of the 
sample are shallowest at the high Peclet numbers. A ratio of the back 
surface gradient to the front surface gradient, which is a measure of the net 
heat loss via solid conduction, is found from Eqn. 1 to be exponentially 
dependent on -Pe. This indicates that for high Peclet number (fast 
regression) burning, despite high heat flux at the burning surface, most of 
that heat is convected back into the flame as the surface regresses, so very 
little heat is actually lost in depth - the ratio is close to zero. On the other 
hand, low Peclet number (slow regression) burning has lower heat flux at the 
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burning surface', but little of that flux is recovered in the' regressive' burning 
and so the ratio of heat flux lost from the hack of the sample to that 
penetrating the burning surface exponentiallv approaches unitv a< Perl<n 
number goes to zero. 

Influences of the Bubble Layer on Temperature Measurements 

The molten polymer and surface bubble layers found in tiit 
experiments were observed to strongly affect measured temperature profiles 
in the solid. Solid phase thermocouples frequently showed strong distortions 
as the regressing surface of the burning PMMA material approached them. 
Some signals dropped precipitously to the glass transition temperature, 
indicating that the thermocouple wire stresses and/or surface tension was 
displacing the bead through the molten viscous polymer to the glass 
transition surface deeper in the burning material (T g = 380 K). This is shown 
clearly in trace A of Figure 5. Within a few mm of the surface, in the highly 
dense bubble region, frequent abrupt changes in the temperature are noted, 
as shown in the B trace. Just above the surface, Trace C shows the gas- 
phase noise associated with the strong vapor jetting at the surface. 



In Fig. !j, two linos are drawn through the average solid and gas-phase 
temperature data. The intersection of these two gradient lines is taken as 
the "surface" temperature (T - (530 Kv although the actual signal from rh<- 
thermocouples very near this location can he quite different, as shown, li i> 
thus very difficult to reproducible and accurately measure local temperature 
gradients from the local thermocouple data due to the influences of the 
molten PMMA and bubble layer. It is frequently necessary to extrapolate the 
surface temperature gradients (gas and solid; from data further removed 
from the surface, which can introduce errors in the calculations. 


Sample Swelling and Regression 

The position of the surface was measured as a function of time during 
each experiment, in order to measure the actual regression rate as a function 
of time. These measurements clearly demonstrate for the first time that the 
sample surface swelled outward significantly during the early phase of the 
experiment. This was attributed to the development of the bubble layer 
during the ignition and solid-phase heatup of the sample. Swelling stopped 
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as the bubble layer became developed and regression began to dominate the 
surface motion. 


An example of this swelling and regression is shown in Figure G for a 
50 cm radius of curvature sample. Ignition time in this experiment is 3 
minutes. During the first few minutes of the experiment after ignition is 
complete, the surface swells outward 4 mm before regression is observable. 
During the swelling period, the flame is gradually moving away from the 
surface as the solid becomes heated and the required conductive flux from the 
flame becomes smaller. Time scales of this period on the order of on the 
solid-phase heatup (~a/V rep . 2 ). As the solid-phase heated layer evolves, so 
does the bubble layer. 

The initial swelling of the sample during and following ignition was 
unexpected. Although other researchers have described a two-phase layer 
during combustion 110 ' 1 '• 21 ", no one has described this initial swelling. 

Thermal expansion of the PMMA as it is heated (estimated to be -0.2 mm) 
is not sufficient to explain the swelling. Rather, the observed swelling 
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appears to he linked to the bubble layer, as is shown in Figure 7, 


The maximum observed swelling in each test is normalized by the 
measured bubble layer thickness and plotted in Figure 7 as a function of 
stretch rate (or sample radius). The surface regression going on during the 
swelling phase is considered to be unrealized swelling since the regression of 
the surface during swelling reduces the measured swelling. The unrealized 
swelling is estimated as shown in Figure G as a linear extrapolation of the 
regression line during the swelling transient. The combined swelling and 
regression normalized by the bubble layer thickness approaches unity as 
stretch rate is reduced, until unstable flame coverage occurs below 3 sec' 1 . 
(For unstable flame coverage, the swelling and regression were not uniform 
across the sample, so these data are not included in Fig. 7.) The increase in 
swelling is thus directly commensurate with the increase in bubbles found in 
these samples during post-burn inspections. 

For larger stretch rates (> 1 2 s' 1 ), the trend in the data indicates that 
small swelling would be masked by regression. This explains why the 
swelling due to the bubble layer has not been observed in previous 
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investigations, where typical buoyant stretch is ai least 20 s '. The burning 
rates in normal gravity are sufficiently high that surface regression is fast 
enough to mask the small amount of swelling that occurs. In addition, tlv 
swelling is so small at the high burning rates due to the thinness of the 
glassy and bubble layers, as shown in Fig. 4. 

Vapor Bubble Motion 

It is obvious that many vapor bubbles reach the surface during the 
experiment. The characteristic sizzling sound of PMMA burning is due to 
these bubbles popping at the surface. In addition, the swelling stops after 
the initial swelling transient, and steady regression is observed. This means 
that there isn't a continual accumulation of the bubbles, so the net 
population of bubbles probably stabilizes during the steady regression 
period. Lastly, there are bubbles left in the sample after it is quenched, so 
the bubbles move so slowly (in either direction) that they do not reach either 
the liquid- gas or solid-liquid boundaries during the quench time (time to 
reach T ? = 105 °C ~ 5 minutes). 
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It would appear from those physical observations that the early period 
of burning has an accumulation of bubbles (swelling) due to the 
ignition/heatup transient, but that the bubble layer later reaches a steadv 
state so steady regression can be observed. This implies that the bubbles are 
formed in depth at approximately the same rate as they reach the surface. 

Buoyancy, thermocapillary convection, and surface regression can all 
cause the bubbles to move relative to the burning surface, as shown in Figure 
8. These three motions define the actual motion of the bubble, as shown in 
Figure 8. The role of buoyancy in the development of the bubble layer is 
worthy of consideration, since the hot gas bubbles will tend to rise through 
the viscous molten plastic. In this experimental configuration, buoyancy 
causes the bubbles to migrate deeper into the solid, away from the hot 
burning surface. The buoyant velocity of a bubble can be estimated to be V b 
~ 1/3 pgR'/p [22] , where p and p are the density and viscosity of the molten 
polymer, R is the bubble radius, and g is the acceleration of gravity. 
Thermocapillary convection causes the bubbles to move toward the hot 
surface along the temperature gradient within the viscous molten material. 
The thermocapillary convective motion of the bubbles can be estimated to be 
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V„ : -Po/PT R VT/3n l ' J| , whore o is surface tension of the* molten polymer and 
T is temperature. Surface regression (V,.,,, r ), measured in each experiment, 
can he considered as a linear coordinate transformation, where a stationarv 
bubble relative to the moving surface can be converted to a moving bubble 
relative to a stationary surface. The regressing surface acts then to move tin* 
bubble toward the burning surface. 

Unfortunately, the viscosity of polymers is a non-iinear function u: 
temperature and molecular weight, and varies from thousands to 10 13 poises 
over the ranges of interest ( "' 31 Similarly, surface tension properties are not 
well characterized. It is thus currently impossible to make a quantitative 
assessment of the bubble motion during the experiment. 

However, recall that the elongated bubbles were aligned normal to the 
curved cylinder's surface rather than along the gravity vector. This 
qualitatively demonstrates the importance of thermocapillary convection 
and/or surface regression on the bubble elongation parallel to the 
temperature gradient within the viscous molten material. If buoyant flow 
were dominant, the elongated bubbles would have all been parallel to each 
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other, and in the direction of the gravity vector, not normal to the curved 
surface. Thus it appears, at least qualitatively, that thermocapillary 
convection and/or surface regression dominate the relative motion between 
the bubbles and the burning surface, and buoyancy plays a minor role. 
Therefore, it is not surprising that these bubble layers have also been found 
in low gravity for comparable low stretch conditions 1 

The glass transition temperature of PMMA is luo C, and the monuiiici 
boiling point is 101°C (1 atm). These two properties dictate that the 
monomer present in the viscous melt will be in the vapor phase. Since 
bubbles appear to generally migrate toward the burning surface, the 
observed bubble layer depth is a conservative measure of the depth of the 
pyrolysis region in the sample. This thickness (3- 7 mm) is quite large when 
compared typical material thicknesses or thermal depths, it is a significant 
fraction of the total effective thickness of the material. 

A common model assumption is that all gasification occurs only at the 
surface via zeroth order Arrhenius kinetics. As was noted earlier, the polymer 
surface is not saturated with monomer, and the monomer formation and 
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"diffusion" to the surface are believed to he the rate-controlling processes of 
gasification 11 ’ 1 . The monomer vapor bubbles' motion to the surface is likely 
to be a significant source of monomer at the surface. Thus an 
understanding of the pyrolysis zone structure of a material requires 
inclusion of the bubble layer characteristics. 


Heatup and Thinning Transients 

Because the flame is very sensitive to the near-surface solid 
temperature gradient, the bubble layer can influence the flame stability. In 
the experiments, the flame stability changes over solid-phase time scales 
associated with sample heatup and surface regression- induced sample 
thinning. During the heatup phase, the sample temperature profile develop 
and the thermal layer develops. For a very near-quenching low stretch case, 
the flame strengthens from flamelets to a uniform flame with gradually 
increasing standoff distance (at 40 min) to a one-dimensional flame with a 
pseudo-stable standoff distance (at 60 min). It remains stable only briefly 
before reverting back into flamelets as the thinning phase begins to 



dominate. The thinning phase is where the thinning of the sample causes 
the temperature gradients to become steeper. These experiment phases arc 
discussed in detail below. 

The experimental solid-phase temperature gradient is shown in 
Figure 9 for the limiting Pe=0.22 run, which was sufficiently weak that a 
one-dimensional uniform flame could be obtained only after the extensive 
transition period from flamelets to a one-dimensional flame. The bubble 
layer/surface temperature gradient measurements agree with predicted 
gradients (described below) for early times (<60 s), but deviate after that. Foi 
Peclet numbers lower than this, a one-dimensional flame was never obtained 
The experimental local temperature gradient at the burning surface for this 
figure is evaluated using a second order polynomial fit to sub-surface 
thermocouples (which effectively move toward the surface in time as the 
surface regresses to expose them). For the low Peclet number of this test, 
the trends in the experimental measurements should be very similar to 
theory if the bubble layer has little effect. 

The transient numerical problem of the conduction-convection energy 
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balance for an unsteady finite thickness material !ls| is solved to gain insight 
into the time scales observed in the experiment. The simple formulation here 
neglects in-depth degradation and property variations due to the presence ->(' 
bubbles, both of which will modify the temperature gradients. Surface 
regression (where L decreases with time), is included, and the predicted 
transient temperature gradients at the burning surface, to which the flame 
is immediately sensitive are compared with measured temperature gradients 
as a function of time, as shown in Figure 9. Conduction into the interior 
decays with time during the heatup transient, but conduction increase again 
at later times as sample thickness decreases due to surface regression. The 
two dominant solid phase transients (heatup and thinning) are such that 
the solid-phase never reaches a steady-state (i.e. never achieves constant 
surface temperature gradient) for a finite-thickness regressing solid. 

The experimental data agree well with the predicted heatup trend, 
showing the heatup transient phasing into the thinning transient. During 
the thinning phase, the local burning surface gradient deviates from the 
anticipated trend at approximately 50 minutes, possibly due to bubble layer 
effects on heat transfer. It is also possible, due to the non-uniform nature of 
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the bubble layer, that the thermocouple point measurements in the bubble 
layer are not representative of the gradient . However, since no anomalous 
thermocouple behavior such as Fig. 5 was apparent, it i< more likelv <ha' it' 
bubble layer development is altering the local gradients. 

The deviation from the predicted trend occurs in conjunction with the 
onset of gas-phase flame sooting also noted at approximately 50 minutes. 

The bubble layer development/swelling continued in this te.->t until at ieast ow 
minutes, at which time surface regression could finally be observed. The 
soot radiation would be expected to significantly increase the overall heat 
flux (25%) , with in-depth absorption of the broad-band soot radiation and 
associated endothermic degradation. The increase in absorbed radiation can 
modify the bubble layer and therefore the local thermal gradients within 
that layer. The measured radiant flux from the surface plus flame, shown in 
the Fig. 9, peaks at approximately 60-65 minutes and is 1.2 W/cm 2 , while 
the surface radiative loss is estimated to be 0.77 W/cm 2 for an emissivity of 
0.85. For an optically-thin gas, the net gas radiation influx thus peaks at 0.43 
W/cm 2 . The total radiant flux from the predominantly blue flames at the 
flamelet transition boundaries (40 min and 85 min) are almost identical to 
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each other at 0.95 W/cnr, or a net of 0.18 W/ctrr from gas-phase radiation. 
The increased gas-phase radiative feedback due to sooting is 0.25 W/cni". 

For comparison, using X. = 2.09 x 10 ! W/cm K the local in-depth eonduetiv<- 
flux varies from 0.7- 1.0 W/cm* during this same period. 

The radiation and flame uniformity are out of phase with the 
predicted solid-phase gradient trend as shown in Figure 9. While the 
minimum gradient occurs at approximately 35 minutes, and increases aiu-i 
that as surface regression (sample thickness) begins to influence the 
gradient. However, the measured radiation peaks later, at 65 minutes, when 
a one-dimensional flame is finally stabilized and the swelling stops. The 
one-dimensional flame is assumed to be stronger than the flamelets that 
preceded it, based upon the radiometer readings shown in Fig. 9 that 
indicate that the radiant emissions are higher for the one-dimensional flame 
than for the flamelets. 

This is counter to what Yang and Then's 1 ' 41 steady theory predicts, 
which is that flame strength decreases as in-depth conduction becomes a 
larger fraction of the net heat transfer from the flame. At 60 minutes, the 
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predicted in-depth conduction has increased , hut the measured gradient is 
hovering near the minimum. If t he local temperat ure gradients at the 
surface are stabilized for an extended period of time due tn the continued 
development of the bubble layer (as indicated by the local gradient data in 
Fig. 9), the gas-phase flame behavior can he explained. As the sample 
regresses further, in-depth conductive loss increases such that the flame 
weakens and stops sooting, which reduces radiative feedback to the solid 
(after 85 min), whereupon the flame weakens back into flamelm.-. Seine 
modeling work addressing the solid phase has already been done [u 21 - v \ but 
a more complete model of flame and solid phase is needed. 

Conclusions 

Large- scale buoyant low stretch stagnation point diffusion flames 
over a solid fuel (polymethylmethacrylate) were studied for a range of 
aerodynamic stretch rates of 2-12 sec' 1 which are of the same order as 
spacecraft ventilation-induced stretch in a microgravity environment. 

An extensive layer of polymer material above the glass transition 
temperature is observed. A bubble layer within the glass layer is found to 
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develop during combustion, and the sample swelling during the heatup 
phase of PMMA combustion is linked directly to the developing bubble layer. 

This swelling, which has not been reported elsewhere i- oplv ea-ilv 
observable when solid-phase Peclet numbers (regression rates) are low. 

The near-surface molten polymer and bubble layers in burning 
PMMA are very complicated. There are property variations with the glass 
transition. The in-depth degradation chemistry processes involved in buobie 
layer development are not fully understood, but in-depth radiant absorption 
and bubble migration may play important roles in the combustion process 
at these low Peclet numbers. 

Thermocapillary convection and surface regression are inferred to be 
the dominant factors in the relative motion between the bubbles and the 
burning surface. Buoyancy appears to play a small-to-insignificant role in 
the bubble migration, so bubble layers are expected to be important in low 
gravity (low stretch) burning of materials. 

Bubble layers in low gravity have previously been shown to be a source 
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of flow (fuol vapor jotting), which enhances burning in a quiescent 
microgravity environment through entrainment of fresh oxidizer. Popping of 
the vapor bubbles at the surface of the burning mat erial has also expelled 
burning globules of the molten material, which poses a fire propagation 
hazard. The bubble layer acted to insulate the polymer surface, which 
stabilized the flame for longer than expected from theory neglecting the 

Despite the ability of these simple model.? to account fur overall soiiu- 
phase behavior, to predict the near-surface bubble layer influences on the 
heat and mass transfer process critical to the dynamic flame stability in 
these low stretch environments, a much more detailed model is needed. 
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Table 1 

PMMA Material Properties 


Property 

Value 

Polymer characteristics 

cast, amorphous, high molecular 
weight 

* i 

Glass Transition Temperature, T„ 

105 C | 

Degradation 

unzips to monomer - C5H802 : 

degradation starts at 200"C 

Monomer Boiling Point 

101'C 
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Table 2 

Bubble Layer Characteristics 


Radius, cm 

5 (half- 
cylinder) 

20 

50 

100 

Stretch Rate, s' 1 

11.7 

5.9 

3.0 . 

2.5 

Ham.- U-\, 

Peclet Number (based on 
initial thickness) 

3.0 

0.55 

0.35 

0.15 

Regression rate, mm/s 
(average) 

0.0087 

0.0029 

0.00125 

0.0009 

flamelets 

bubble layer thickness, 
mm 

3.5 

4.5 

5 

7 

bubble diameter, typical, 

mm 

.2 - 2 

0.4-1 

0.5- 1.5 

i 

u.5-o 

bubble layer character, in 
brief 

larger, pseudo- 
spherical 
bubbles deeper 
in layer, many 
smaller at surf. 

very 
dense, 
elongate 
d glassy 
bubbles 

large 

crinkled 

opaque 

tethered 

bubbles 

larger 

crinkled 

opaque 

tethered 

bubbles 

surface crater sizes, 
typical, mm 

0.3 

0.5-2 

0.5-5 

0.5-5 

surface 'color' (exclusive of 
black drip areas) 

black 

light 

caramel 

caramel 

caramel 

refractive index surface 
relative to front surface, 
mm 

10 

12 

9.5 

11 

final sample 
thickness, mm 

47 

19 

13 

17 

swelling observed, mm 

1 

3.5 

4 

5 

time to peak swelling, s 

100 

400 

600 

600 

back surface cooling 

ambient 

ice 

ice 

ice 
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I iguro 1 . Schematic concept oi the edge view ol a burning cylindrical PMMA sample. 

l igiire 2: Magnihed images ol the re-soliditied burned surface structure and opacity for different 
radii PMMA samples: a)5 cm, b)20 cm, c) 50 cm. d) 100 cm. 

figure 3: Magnified images of the cross section view of the re-solidified PMMA samples 
showing the bubble layers for different radii: a)5 cm. b)20 cm. c) 50 cm. d) 100 cm. 

figure 4: Steady-state non-dimensional temperature profile solutions for various Peelet 
numbers. Y=0 at the burning surface. Experimentally measured bubble layer thickness (•) 
and the position of the index of refraction surfaco(A) compare favorably with theorv. 
0g is the non-dimensional glass transition temperature. Data for ice bath back surface cooline 
tests only. 

Figure 5: Examples of solid-phase thermocouple traces showing anomalous behavior as the 
regressing PMMA surface approaches the thermocouple bead. 

Figure 6: Time history ot the location ot the burning PMMA surface during a test, showing 
swelling and regression behavior for Pe=0.35 case (50 cm). Ignition time is 3 minutes, after 
which a stable flame is maintained over the sample. 

Figure 7: Swelling, surface regression distance during swelling, and combined swelling and 
regression, normalized by the measured bubble layer thickness (right axis) as a function of the 
stretch rate. The combined value approaches unity at low stretch, indicating the bubble laver 
thickness is directly responsible for the observed swelling. Data for stable one-dimensional 
flames only. 

Figure 8: Schematic of the forces acting on a bubble within the molten polymer layer. 

Figure 9: A comparison of experimental and theoretical solid-phase temperature gradients alona 
with radiant flux as a function of time for Pe=0.22 case (75 cm), showing good agreement during 
the heat up phase, but departures during the thinning phase that are attributed to the bubble layer. 
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